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The spectral properties of hybrid meson interpolating fields are investigated. The quantum numbers of the 
meson are carried by smeared-source fermion operators and highly-improved chromo-electric and -magnetic field 
operators composed with APE-smeared links. The effective masses of standard and hybrid operators indicate 
that the ground state meson is effectively isolated using both standard and hybrid interpolating fields. Focus is 
placed on interpolating fields in which the large spinor components of the quark and antiquark fields are merged. 
In particular, the effective mass of the exotic 1 _+ meson is reported. Further, we report some values for excited 
mesonic states using a variational process. 



1. INTRODUCTION 

Major experimental efforts are currently aimed 
at determining the possible existence of exotic 
mesons; mesons having quantum numbers that 
cannot be carried by the minimal Fock space com- 
ponent of a quark-antiquark pair. Of particular 
mention is the proposed program of the GlueX 
collaboration associated with the forthcoming up- 
grade of the Jefferson Laboratory facility. The 
observation of exotic states and the determina- 
tion of their properties would elucidate aspects of 
QCD which are relatively unexplored. 

The quantum numbers J PC = H , 0~~, 1 
etc. cannot be carried by a quark-antiquark pair 
in a ground-state S'-wave. Lattice QCD calcu- 
lations exploring the non-trivial role of explicit 
gluon degrees of freedom in carrying the quan- 
tum numbers of the meson suggest that exotic 
meson states do indeed exist and have a mass the 
order of 2 GeV £Q. These findings are further 
supported here. 

2. SIMULATION METHODOLOGY 

Operators carrying exotic quantum numbers 
can be constructed by merging standard local 
interpolating fields q a (x)Tq a (x) with chromo- 
electric, Ef b (x), or chromo-magnetic fields, 



Bf h (x). The J PC quantum numbers of the in- 
terpolator are derived from the direct product 
of those associated with the quark bilinear and 
Ef {1--) or Bf (l + ~). For example, com- 
bining the vector current of the p meson with 
a chromo-magnetic field, 1 ® l" 1 provides 
0"+ © 1-+ © 2"+ with the 0"+: q a j l q b Bf (it 
meson) and the 1 h eijkq a liq h B°' b (exotic). We 
restrict ourselves to the lowest energy-dimension 
operators, as these provide better signal with 
smaller statistical errors. Tabled summarizes the 
standard and hybrid interpolating fields explored 
herein. 

The formulation of effective interpolating fields 
for the creation and annihilation of exotic meson 
states continues to be an active area of research. 

For example, one can generalize the structure of 
the interpolating fields further to include nonlocal 
components where link paths are incorporated to 
maintain gauge invariance and carry the nontriv- 
ial quantum numbers of the gluon fields. In this 
case, numerous quark propagators are required 
for each gauge field configuration rendering the 
approach computationally expensive. 

Here we consider local interpolating fields. 
Gauge-invariant Gaussian smearing |2I3| is ap- 
plied at the fermion source (f = 3), and lo- 
cal sinks are used to maintain strong signal in 
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Table 1 



J quantum numbers and their associated meson interpolating fields. 



0++ 


o+- 


o-+ 


o— 


-a q a 

-iq a jjEfq b 
-q a ljl4l 5 Bfq b 
-q a ljHEf>q h 


q a liq a 
q a l5l] Bfq b 


q a l5 q a 

q a i5nq a 

~~q a lj Bfq b 
-q a j4ljBfq b 


-iq a l5lj Efq b 




1++ 


1+- 


1-+ 


I — 


-iq a i5ij<f 

iq a ^Bfq b 
ie ]kl q a lkEfq b 
ie jk iq a lkl4E? b q b 


-iq a i5i4ij(f 

iq a Bfq b 
q a lb Efq b 
q a l5li Efq b 


q a j 4 Efq b 

-a„, r>ab „b 
Sjkiq 7fc-D; 9 
r>ab~b 

tjkiq j4jk-&i q 

-iejkiq a l5l4lkEfq b 


—iq ar )jq a 
q a Efq b 
-iq a ~f 5 Bfq b 
iq a lAl5 Bfq b 



the two-point correlation functions. Chromo- 
electric and -magnetic fields are created from 
APE-smeared links [3j at both the source and sink 
using the highly-improved 0(a 4 )-improved lattice 
field strength tensor |S] . In this study, the smear- 
ing fraction a = 0.7 (keeping 0.3 of the original 
link) and the process of smearing and SU (3) link 
projection is iterated four times jSj. This amount 
of smearing is sufficient to provide a meaningful 
topological charge and appears to be suitable for 
the creation of exotic mesons. As such, the results 
presented here supersede an earlier presentation 
of hybrid meson masses . 

Propagators are generated using the fat-link ir- 
relevant clover (FLIC) fermion action 8 where 
the irrelevant Wilson and clover operators of the 
fermion action are constructed using fat links 
while the relevant operators use the untouched 
(thin) gauge links. FLIC fermions provide a new 
form of nonperturbative 0(a) improvement |9I1U| 
where near-continuum results are obtained at fi- 
nite lattice spacing. Access to the light quark 
mass regime is enabled by the improved chiral 
properties of the lattice fermion action |10j . 

Excited states are extracted using a variational 
technique, corresponding to a construction of op- 
timal linear combinations of the original opera- 
tors. For the sake of completeness, we shall dis- 
cuss this here, in direct analogy to the procedure 
described in 



3. ANALYTICAL PROCESS 

Consider the momentum-space meson two- 
point function for t > 0, 

G l3 (t,p) = ^e-* s (n\ Xi {t,g)rfi{oMn) > (!) 

X 

where i,j label the different interpolating fields 
and we focus on Lorentz scalar interpolators for 
simplicity. At the hadronic level, 



x (jT,p'| x t(o,a)|n> , 

where the \H,p') are a complete set of hadronic 
states. 



J2\H,p')(H,p'\=I . 

H,p' 



(2) 



We can make use of translational invariance to 
write this as 



E 



x,H.p f 



n 



H,p'\x](0) 



X*(0)e 



n 



H,p') x 



= J2e- EHt (n\Xi\H,p)(H,p\x]\n 



(3) 



H 



It is convenient in the following discussion to la- 
bel the states which have the x interpolating field 
quantum numbers as \H a ) for a = 1, 2, • • • , N. In 



3 



general the number of states, N, in this tower of 
excited states may be infinite, but we will only 
ever need to consider a finite set of the lowest 
such states here. After selecting zero momentum, 
p = 0, 

N 

G iS (t) = Gij (i, 0) = £ e-^'Af A}° , (4) 

a=l 

where Xf and X^ a are coefficients denoting the 

couplings of the interpolating fields Xi an d x\i 
respectively, to the state \H a ). If we use identical 
source and sink interpolating fields then it fol- 
lows from the definition of the coupling strength 
that A]" = (Xf)* and from Eq. Q we see that 
Gij(t) = [Gji(i)]*, i.e., G is a Hermitian matrix. 
If, in addition, we use only real coefficients in the 
link products, then G is a real symmetric ma- 
trix. For the correlation matrices that we con- 
struct we have real link coefficients but we use 
smeared sources and point sinks and so in our cal- 
culations G is a real but non-symmetric matrix. 
Since G is a real matrix for the infinite number 
of possible choices of interpolating fields with real 
coefficients, then we can take Xf and X\ a to be 
real coefficients here without loss of generality. In 
constructing correlation functions, we effectively 
average over {U} and {U*} configurations to en- 
sure A" is purely real, even on a finite ensemble 
of gauge field configurations [TT] . 

Now, let us consider the ideal case where we 
have N interpolating fields with the same quan- 
tum numbers, but which give rise to N linearly 
independent states when acting on the vacuum. 
In this case we can construct TV ideal interpolat- 
ing source and sink fields which perfectly isolate 
the N individual hadron states \H a ), i.e., 

N 

0t« = 5>? X J, (5) 
i=i 

N 

r = $>rxi, (6) 

i—l 

such that 

{H p \^ a \n) = S a0 z^ a , (7) 

<t> a \H P ) = 5 a0 z a , (8) 



where z a and z' a are the coupling strengths of <p a 
and (j>' !a to the state \H a ). The coefficients uf and 
v* a in Eqs. J5J may differ when the source and 
sink have different smearing prescriptions, again 
indicated by the differentiation between z a and 
z* a (recall z is real). 

For notational convenience for the remainder 
of this discussion repeated indices k are to 
be understood as being summed over, whereas a 
denoting a particular state is not. At p — 0, it 
follows that, 

= Xfz^e-™"*. (9) 

The i-dependence in this expression is purely in 
the exponential term, leading to the recurrence 
relationship 

Gij(t)uf = e m <*G ik (t + l)u% , (10) 
which can be rewritten as 

[G(t + l)]ltG ij (t)uf = e™«u a k . (II) 

This is the generalized eigenvalue equation for 
[G(t + l)]~ 1 G(t) with eigenvalues e m ° and eigen- 
vectors u a . Hence the natural logarithms of the 
eigenvalues of [G(t + l)]~ 1 G(t) are the masses of 
the N hadrons in the tower of excited states for 
the given quantum numbers . The eigenvectors 
are the coefficients of the \ fields providing the 
optimal linear combination for that state. 

One can also construct the equivalent left- 
eigenvalue equation to recover the v vectors, pro- 
viding the optimal linear combination of annihi- 
lation interpolators, 

v% a G kj (t)=e m «v*<*G ij (t+l) . (12) 
Recalling Eq. JHJ, one finds: 

Gij(t) uf = z^XTe- m ^ , (13) 
v* a Gij(t) = z a Xfe- m « t , (14) 
v* k a G kj (t)G u (t) uf = z a z^ a XfX] a e~ 2m ^(15) 
The definitions of Eqs. Q and |JEJ) imply 
v* a G i:i (t) uf = z a z fa e- m " t , (16) 
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indicating the eigenvectors may be used to con- 
struct a correlation function in which a single 
state mass m^is isolated and which can be an- 
alyzed using the methods of Section II. We refer 
to this as the projected correlation function in the 
following. Combining Eqs. 115fl and 116(1 leads us 
to the result 



vt a G kj (t)G a (t) uf 
vt a G kl (t)uf 



A? A 



(17) 



By extracting all iV 2 such ratios, we can exactly 
recover all of the real couplings Af and \l a of Xi 

and x} respectively to the state \H a ). 

Note that throughout this section no assump- 
tions have been made about the symmetry prop- 
erties of Gij. This is essential due to our use of 
smeared sources and point sinks. 

In practice we will only have a relatively small 
number, M < N, of interpolating fields in any 
given analysis. These M interpolators should be 
chosen to have good overlap with the lowest M 
excited states in the tower and we should attempt 
to study the ratios in Eq. 1(17(1 at early to inter- 
mediate Euclidean times, where the contribution 
of the (N — M) higher mass states will be sup- 
pressed but where there is still sufficient signal to 
allow the lowest M states to be seen. This pro- 
cedure will lead to an estimate for the masses of 
each of the lowest M states in the tower of ex- 
cited states. Of these M predicted masses, the 
highest will in general have the largest system- 
atic error while the lower masses will be the most 
reliably determined. Repeating the analysis with 
varying M and different combinations of interpo- 
lating fields will give an objective measure of the 
reliability of the extraction of these masses. 

In our case of a modest 2x2 correlation matrix 
(M = 2) we take a cautious approach to the selec- 
tion of the eigenvalue analysis time. As already 
explained, we perform the eigenvalue analysis at 
an early to moderate Euclidean time where sta- 
tistical noise is suppressed and yet contributions 
from at least the lowest two mass states is still 
present. One must exercise caution in perform- 
ing the analysis at too early a time, as more than 
the desired M = 2 states may be contributing to 
the 2x2 matrix of correlation functions. 
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Figure 1. Effective mass plot of the 1 + ex- 
otic meson obtained from the hybrid interpolating 
field e 3k iq a lk Bfq b . 



4. RESULTS 

The following results are based on 110 mean- 
field O(a 2 )-improved Luscher-Weisz gauge 
fields on a 16 3 x 32 lattice at j3 = 4.60 provid- 
ing a lattice spacing of a = 0.122(2) fm set by the 
string tension ^fa = 440 MeV. 

Of the hybrid interpolators listed in Table ^ 
only the interpolating fields merging the large 
spinor components of the quark and antiquark 
fields provide a clear mass plateau. The effective 
mass plot for the exotic 1 ^ meson is illustrated 
in Fig. ^ where a plateau at early times is ob- 
served confirming the existence of the exotic 1 h . 

Figures [21 and [21 illustrate the effective masses 
M(t) = -log(G(t + 1)/G(f)), obtained from 
the first and third, and second and fourth, pion 
(0 interpolators of Table 2] for our interme- 
diate quark mass (ml ~ 0.6 GcV 2 ). Excel- 
lent agreement is seen between the standard and 
hybrid interpolator-based correlation functions. 
Similar results are seen in Fig. 0] comparing ef- 
fective masses obtained from the first and third 
p-meson (1~~) interpolators of Tabled 

Table [2 summarizes these preliminary results, 
and table [3| presents a further preliminary calcu- 
lation of a pion excited state mass, using the first 
two pion operators listed above. 

Further work on this topic will focus on increas- 
ing the statistics, and increasing the number of 
operators used in the variational process. 
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Table 2 

Meson masses as a function of the hopping parameter k. 



Mass(GeV) 







Operator 


k = 0.1260 


K = 0.1266 


k = 0.1273 


k = 0.1279 


k = 0.1286 


7T 


o-+ 




0.965 ± 


.006 


0.887 ± 


.006 


0.789 ± 


.007 


0.696 ± 


.007 


0.570 ± 


.008 








0.957 ± 


.005 


0.879 ± 


.005 


0.781 ± 


.005 


0.688 ± 


.006 


0.563 ± 


.006 






-q a f j Bfq b 


0.998 ± 


.030 


0.916 ± 


.031 


0.813 ± 


.031 


0.717± 


.032 


0.589 ± 


.035 






-q a lilj Bfq h 


0.984 ± 


.033 


0.902 ± 


.034 


0.800 ± 


.035 


0.704 ± 


.037 


0.575 ± 


.039 


h 


1+- 




1.713 ± 


.018 


1.671 ± 


.019 


1.623 ± 


.021 


1.583 ± 


.022 


1.541 ± 


.026 






iq a Bfq b 


1.685 ± 


.205 


1.621 ± 


.226 


1.525 ± 


.275 


1.398 ± 


.353 


1.120± 


.500 


P 


1— 


—iq a jjq a 


1.212 ± 


.009 


1.157 ± 


.010 


1.093 ± 


.012 


1.037 ± 


.014 


0.973 ± 


.019 






q a Efq h 


1.198 ± 


.056 


1.151 ± 


.063 


1.099 ± 


.076 


1.048 ± 


.091 


0.958 ± 


.110 






-iq a j 5 Bfq b 


1.214 ± 


.050 


1.156 ± 


.053 


1.085 ± 


.058 


1.018 ± 


.064 


0.922 ± 


.077 






iq a lil5 Bfq b 


1.125 ± 


.049 


1.060 ± 


.054 


0.982 ± 


.061 


0.907 ± 


.069 


0.799 ± 


.087 




1-+ 


-tjkiq a JkBf b q b 


2.519 ± 


.186 


2.483 ± 


.192 


2.450 ± 


.200 


2.434 ± 


.207 


2.430 ± 


.211 



Table 3 

Pion excited state from variational analysis. 





k = 0.1260 


Mass(GeV) 
k = 0.1266 k = 0.1273 


k = 0.1279 


K = 0.1286 


tt(I): 
tt(2): 


0.956 ± .007 
1.889 ± .064 


0.878 ±.007 0.779 ±.008 
1.819 ±.070 1.730 ±.080 


0.684 ± .009 
1.647 ± .096 


0.555 ± .014 
1.548 ± .141 
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Figure 2. Effective mass plot for correlation func- 
tions of the standard pion interpolator q a j5q a and 
the hybrid pion interpolator q a jjBj b q b . 
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Figure 3. Effective mass plot for correlation func- 
tions of the standard axial-vector pion interpo- 
lator <? a 7574(? a and the hybrid pion interpolator 

q a -f4j-jBfq b . 
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